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FOREWORD 


Xhls  report  vat  prepared  by  Or.  Dennis  S.  Berry  under 
the  direction  of  Charles  Falrhurst  of  the  School  of  Mineral 
and  Metallurgical  Engineering,  University  of  Minnesota, 
Minneapolis,  Minnesota  under  Contract  DA-25-066-ENQ  14,765 
with  th<  Missouri  River  Division,  Corps  of  Engineers. 

Dr.  Dennis  S.  Berry  of  the  University  of  Nottingham, 
Nottingham,  England  was  principal  investigator  for  this 
research. 

L.  B.  Underwood,  Division  Geologist  monitored  this 
research  contract  for  the  govern  sent;  K.  S.  Lane  was  Chief 
of  Oeolt;;{y,  Soils  and  Materials  Branch,  and  J.  0.  Ackerman 
was  Chief  of  Engineering  Division  of  the  Missouri  River 
Division,  Corps  of  Engineers  during  this  investigation. 

Funds  were  provided  by  the  Office,  Chief  of  Engineers, 
Department  of  the  Amy  under  Military  Construction  and 
Investigational  Programs,  0AM,  7.60.12A  QB- 1-02-007. 

Based  on  the  findings  of  this  report  continuing 
research  is  being  sponsored  by  Army  0AM  funds. 


PREFACE 


This  report  is  part  of  a  continuing  study  into  Methods  of 
Determining  In-Situ  Rock  Stresses  at  Great  Depths.  Results  of  the 
first  part  of  the  study  were  reported  in  Technical  Report  No.  1-68, 
bearing  the  title  underlined  above  and  published  in  February,  1968 
by  the  Missouri  River  Division,  Corps  of  Engineers,  Omaha.  In  that 
report  it  was  noted  that  all  current  techniques  of  in-situ  rock  stress 
determination  were  based  on  the  assumption  that  the  rock  behaved  as 
a  linearly  elastic,  isotropic,  continuum.  Since  many  rocks  are  not 
isotropic  it  was  felt  that  an  analysis  of  the  influence  of  rock  anisotropy 
on  the  accuracy  of  the  techniques  was  needed.  It  was  therefore  decided 
to  attempt  solutions  to  problems  in  whicn  the  rock  was  considered  to 
behave  as  a  transversely  isotropic  elastic  material.  A  transversely 
isotropic  material  is  one  for  which  the  elastic  properties  are  invariant 
with  respect  to  rotations  about  a  single  axis  only  in  the  material.  Five 
independent  elastic  constants  (see  Technical  Report  1-68  p.  6)  are 
necessary  to  define  a  transversely  isotropic  material  compared  to  two 
for  isotropic  material.  Bedded  or  laminated  rocks  such  as  shales  and 
gneisses  appear  to  be  fairly  well  described  by  the  transversely  isotropic 
model.  The  use  of  more  sophisticated  models,  involving  a  greater  number 
of  physical  constants ,  -was  considered  unwarranted  since  experimental 
determination  of  the  constants  would  be  difficult  and  the  mathematics 
quickly  becomes  intractable. 

The  theoretical  solution  for  the  Stresses  on  the  Surface  of  a 
Circular  Hole  in  an  Infinite  Transversely  Isotropic  Elastic  Medium  due 
to  General  Stresses  at  Infinity  and  Hydrostatic  Pressure  at  the  Hole,  was 
obtained  by  Dr.  D.  S.  Berry,  Department  of  Theoretical  and  Applied 
Mechanics,  University  of  Nottingham,  Nottingham,  England,  and  is 
presented  in  Appendix  1  of  Technical  Report  1-68.  The  solution  enables 
the  influence  of  rock  anisotropy  to  be  considered  in  the  analyst,  of  the 
hydraulic  fracturing  technique  of  stress  determination. 

This  report,  The  Theory  of  Stress  Determination  by  Mear  s  of 
Stress  Relief  Techniques  in  a  Transversely  Isotropic  Medium  ,  also  by 
Dr.  Berry,  presents  the  analytical  expressions  for  the  strain.'  and  dis¬ 
placements  at  the  surface  of  a  circular  hole  in  an  infinite  tr;  nsversely 
isotropic  elastic  medium  due  to  general  stresses  at  infinity,  This  solution 
permits  the  influence  of  rock  anisotropy  to  be  considered  i  i  the  stress- 
relief  (overcoring)  tecnniques. 

Investigation  of  the  effect  of  rock  anisotropy  on  stresses  in  an 
elastic  inclusion,  the  theoretical  basis  of  tne  remaining  mportant  class 
of  stress-determination  methods,  is  now  in  progress  and  will  be  published 
in  a  subsequent  report. 

C.  Fairhurst 

September  6,  1968 


THE  THEORY  OF  STRESS  DETERMINATION  BY  MEANS  OF 
STRESS  RELIEF  TECHNIQUES  IN  A  TRANSVERSELY 
ISOTROPIC  MEDIUM 

D.  S.  Berry 


1 .  Introduction 


The  use  of  the  stress  relief  technique  requires  theoretical 
knowledge  of  the  radial  displacement  in  a  long  hole  as  the  stress  is 
removed  by  overcoring  or  other  means.  Panek  (1966)  has  pointed  out 
that  relief  of  stress  in  the  axial  direction  affects  the  radial  displace¬ 
ment  in  an  isotropic  medium,  while  Berry  and  Fairhurst  (1966)  have 
incorporated  this  effect  in  calculating  results  for  a  transversely 
isotropic  medium  in  which  a  hole  is  drilled  in  the  plane  of  symmetry 
or  normal  to  it.  This  report  extends  those  results  to  an  arbitrarily 
oriented  hole  drilled  into  a  medium  with  any  homogeneous  state  of 
stress . 


The  author's  previous  work  on  the  stress  around  a  hole  in  a 
transversely  isotropic  medium  (Appendix  I  of  Tech.  Rpt.  No.  1-68 
by  C .  Fairhurst,  referred  to  here  as  I)  is  drawn  upon  freely  and  the 
same  notation  is  adopted.  The  method  of  solution  is  based  upon  the 
work  of  Milne-Thomson  (1962)  on  "anti-plane"  strain. 

2.  Displacements  due  to  antiplane  strain 

As  in  I,  the  plane  z  =  0  defines  an  arbitrary  cross-section 
of  the  hole,  the  x-axis  is  chosen  to  lie  in  a  plane  of  elastic  symmetry 
(the  lines  of  intersection  of  these  planes  with  the  plane  of  the  paper 
are  indicated  in  Figs.  1  and  3  )  in  such  a  way  that  the  angle  0  that 
the  plane  makes  with  the  y  -  axis  satisfies  the  condition  0  5  0  5  tt/2 
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while  a  is  obtained  from  equation  1(2 . 5)  : 
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6,  =  +  k0-3a  +  koo<*  +  ko/.  T  . 

z  13  x  23  y  33  z  34  yz 


(2.2) 


The  constants  n  ,  ,  n3  are  given  by  1(2.20). 

Substitution  of  the  expressions  (2.1)  and  (2.2)  in  the  relations 
1(2.6)  (valid  for  constant  %  )  gives  the  following  expressions  for  the 
strain  components: 
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and  the  -trg  and  are  given  by  I  (2.6  ). 

If  we  write  V  (z  )  =  J  W  (z  )  dz  *  then  integration 
of  (2.3)  with  respect vto^x  and  (2.v4)  vwithv respect  to  y  yields 
the  following  two  expressions  for  components  of  displacement 
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where  f  (y)  and  f„(x)  are  arbitrary  functions.  An  expression 
for  2y  =  du/dy^+  dv/dx  can  be  found  from  (2.9)  and  (2.10) 
and  comparison  with  (2.5)  shows  that  fj(y)  and  f2(x)  are  constants 
(which  may  be  taken  as  zero  since  they  correspond  to  rigid  body  motion 
only) . 

Since  we  are  considering  deformations  which  are  independent  of 
z,  we  have  that 
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Then  integration  of  (2 . 6)  and  (2 . 7)  give 
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Hence ,  apart  from  constants ,  which  may  be  ignored 


f3W  s  0  ,  f4  (y)  =  K4ezy  . 


Also,  we  find  that 
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where  d  (X)  and  d  (X)  are  given  by  1(2.18)  and  this  is  identically 
zero  by  the  definition  of  in  1(2.20). 

Collecting  the  above  results  we  have  the  following  expressions 
for  the  three  components  of  displacement: 
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Formulation  of  the  problem. 


In  the  stress-relief  method,  diametral  displacement  indicators 
(possibly  strain-gauges  also)  are  inserted  at  a  section  of  a  hole 
sufficiently  remote  from  its  ends  for  their  influence  to  be  negligible, 
the  surrounding  stress  is  relieved  and  changes  in  diametral  measure¬ 
ments  of  the  hole  are  recorded,  possibly  changes  in  axial  strain  also. 

If  the  regional  stress  components  are  denoted  by  (  oij)0  ,  as  in  Section  3 
of  I,  we  can  describe  their  removal  as  the  application  of  stress  com¬ 
ponents  -  (  0^j)o  #  and  the  boundary  conditions  can  be  stated  as 
follows: 


°r  *  Tre  =  Trz  =  0  '  onr*a  ' 


(3.1) 


a  -  -  (a  )  at  °° 
ij  v  ij  *o 


Now  these  conditions  are  just  those  of  1(3.1)  apart  from  the  minus  sign 
in  front  of  (  Ojj  )Q  .  In  addition,  it  is  specified  in  I  Section  3  that 
ez  should  be  constant  everywhere:  if  we  can  show  that  that  condition 
is  satisfied  in  the  stress-relief  problem  then  the  stress  functions  for 
the  solution  of  1(3.1)  can  be  adopted  for  the  present  problem,  merely 
with  a  change  of  sign.  It  can  be  demonstrated  as  follows. 


Since  we  are  assuming  that  the  material  conforms  to  the 
assumptions  of  classical  elasticity  we  can  invoke  the  principle  of 
superposition  and,  in  particular,  we  know  that  the  chronological 
order  in  which  two  or  more  constraints  are  imposed  does  not  affect 
the  resulting  elastic  field.  Consequently,  the  axial  strain  in  a  stress- 
relieved  sample  is  the  same  as  if  the  hole  were  made  after  removal  of 
the  sample  instead  of  before.  Now,  at  a  cross  section  remote  from  the 
ends  we  can  safely  assume  that  the  in  situ  drilling  produces  no  axial 
strain.  Also,  drilling  a  hole  after  removal  produces  no  axial  strain 
because  the  sample  now  has  no  applied  stresses  capable  of  producing 
any  deformation.  However,  the  final  state  is  the  same  in  each  case, 
so  the  axial  strain  is  just  that  which  is  produced  by  relieving  stress 
in  the  sample  without  a  hole.  Thus  the  axial  strain  is  not  dependent 
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on  the  presence  of  the  hole  and  so  is  independent  of  position  in 
the  rock. 

If  we  write  out  the  second  set  of  conditions  in  (3. 1)  in  non¬ 
tensor  notation,  using  dj,  ....  t^,  ...  for  the  regional  stress 
components,  we  have  at  » 
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If  these  values  are  substituted  in  (2.2)  we  obtain  the  axial  strain 

in  terms  of  the  regional  stress  components,  valid  for  all  x  and  y 

since  e  is  constant: 
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which  is  just  1(3.2),  apart  from  the  change  of  sign  giving 
etc. 
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4.  Solution  of  the  general  problem. 

The  conditions  to  be  satisfied  are  the  values  of  the  boundary 
tractions  and  stress  components  at  infinity,  given  by  (3.1),  together 
with  the  general  condition  on  the  stress  components  given  by  (3.3). 

The  complex  potential  method  gives  stress  solutions  in  the 
form  (2.1),  apart  from  the  component  az  which  can  be  obtained  from 
relation  (3.3). 

The  solution  of  1(3.1)  under  condition  1(3.2)  gives  stress 
components  denoted  by  (a^)j  etc .  in  I  and  is  obtained  from  the 
complex  potentials  denoted  by  [W  (z  )  ]1#  Clearly  the  solution  given 
by  complex  potentials 
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will  give  stress  components  everywhere  and  o  =  -(or  )q 

at  infinity  in  particular,  thus  satisfying  the  conditions  (3.1)  and  (3,3). 

Thus,  from  Section  4  of  I  the  stress-relief  solution  is  given  by 
(3.2)  and  (3.3)  and  the  expressions 
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The  various  constants  in  (4.2)  to  (4.5)  are  given  by  equations  1(4.6)  to 
I  (4.9),  with  the  addition 
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with  the  A^  given  by  1(4. 12)  ,  while  it  is  easy  to  show  that 
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Substitution  of  (4.6)  and  (4.7)  into  (4.3)  leads  to  the  expression 
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and  so  by  integration  we  find  that 
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For  the  stress-relief  problem  we  require  the  radial  displacement 
at  the  boundary  of  the  hole  ,  r  =  a .  At  the  boundary 
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(The  choice  of  sign  here  is  immaterial  because  of  the  symmetry  of  the 
solution) . 

Substitution  in  (4. 10)  leads  to  the  result  for  r  =  a 
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To  obtain  the  radial  displacement  from  these  expressions  and 
the  displacement  expressions  (2.13)  we  note  first  that 
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r0  /  .  lvX  vovv 

v=  1  v 


4  (L  4  1  ~~  )  ( a  Y  e"219-  A  )  ] 
lv  -  VO  V  v  1 

V 


4  eze  iS  (K1cos04iK2sin9  )'J. 


We  now  put 


L.  L, 

G  =  L.  +  i  — ^  ,  H  =  L.  -  i 
v  lv  X  v  lv  X 

v  v 


and  take  the  real  part  of  (4.1 5)  to  obtain 


u 


=  A  4  B^cos  29  C  sin29  , 
a  0  0  0 


where 


A  =  Re  L  (AH  *  a  Y  G  )  |  e  (K  +  KJ  . 

0  1  V  V  vovv  *  z  1  2 


y 


#  Re  (AG  -  a  y  H  )  -  £e  (K  -  O 
0  ,  v  v  vo  v  v  *  z  1  2 


v-  1 
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L 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


Im  Y  (A, 

v=  1 


G  +  a  Y  H  ) 

V  V  VO  V  V 


Equations  (4.17)  and  (4.18)  give  tne  diametral  expansion  of 
the  nole  (wnen  tne  regional  stress  is  relieved)  in  terms  of  tne  regional 
stress  components  (through  tne  jv0  and  ez  )  and  the  elastic  constants 
of  the  medium.  Tne  orientation  of  tne  hole,  0  ,  determines  the  values 
of  the  constants  wnich  actually  appear  in  tne  above  equations,  so  that 
the  values  of  A^  ,  B^and  are  dependent  on  0  as  well  as  on  the 
regional  stress  components  for  a  given  rock  mass. 

5  .  The  exceptional  cases. 

As  the  author  explained  in  I,  the  cases  0  -  0,  n/2  are  not 
covered  by  tne  general  solution  for  antiplane  strain  in  a  transversely 
isotropic  medium  because  the  stress  functions  \  (x,y)  and  |  (x,y) 
then  satisfy  independent  equations,  instead  of  the  mixed  equations 
1(2.9)  and  1(2.10).  The  result  is  that  components  of  displacement  u 
and  v  (or  ur  and  ug  )  are  independent  of  the  regional  stress  components 

T13  anci  T23  ’ 

The  required  results  for  radial  deformation  due  to  stress  relief 
have  been  given  by  Berry  and  Fairhurst  (1960)  in  terms  of  principal 
stress  components  assumed  to  lie  parallel  to  the  axis  of  the  hole  and 
normal  to  it.  The  derivation  of  these  results  is  given  in  this  section, 
using  the  notation  of  this  report. 

The  result  for  an  isotropic  material  can  be  easily  obtained  from 
either  of  these  results. 

5 . 1  Case  0  =  0 

We  first  show  the  close  correspondence  between  plane- 
strain  problems  in  this  case  and  similar  problems  in  an  isotropic 
material.  This  enables  us  to  use  the  well  known  plane-strain  result 
for  a  stress-free  hole  in  an  isotropic  material.  The  correction  for 
non-zero  axial  strain  is  then  applied  to  complete  the  solution. 

From  1(1 . 5)  and  1(5 . 1) ,  the  stress-strain  relations  for  0  0 

reduce  to 
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e  =  b,  d  •  b.0o  -f  b  o  , 
x  1 1  x  12  y  13  z 


e  ^  b,  „o  +  b,  ,a  b,  .a 
y  12  x  11  y  13  z 


e  -  b  „a  -  b  a  -  b0.o  , 
z  13  x  13  y  33  z 


2y  =  b  .t 
yz  44  yz 


2y  -  bAAT 
zx  44  zx 


2y  =  2  (b  -  b  J  t 
xy  11  12  xy 


where  the  coefficients  are  given  by  1(1.3)  in  terms  of  moduli  and 
Poisson's  ratios . 

For  plane -strain,  e  =  y  =  y  =0  and  hence 
*  z  yz  zx 


(5.1) 


a  =  -  (  a  a  ) . 

2  b33  x  y 


(5.2) 


By  substitution  from  (5.2)  the  remaining  three  equations  then  become 


€  .0  '  .  i*0  , 

x  11  x  12  y 


e  =  t,  0  a  -l  a  , 
y  12  x  11  y 


y  =  {t.  -  l.  J  T 

xy  11  12  xy 


(5.3) 


where 
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JStJ 


'•*?*"■**»'«*-  r 


12 


E2V2 


'1 


(5.4) 


using  1(1 . 3) . 


Now  equations  (5.4)  nave  precisely  the  same  form  as  the 
plane-strain  equations  of  an  isotropic  medium.  It  follows  that  the 
same  methods  of  solution  are  applicable  (we  have  already  seen  in  I 
that  the  stress  function  is  biharmonic),  and  the  actual  solutions  are 
identical  so  long  as  the  constants  of  (5.4)  are  substituted  for  the 
corresponding  constants  for  an  isotropic  material.  The  latter  are 


,1  =  1  ~  v 

'll  E 


1  -  v  ,  I 

~27  and  1  12 


Solving  for  p  and  v  we  have 


b  = 


1 


2(tn-5I2) 


and  v  -  - 


v(l  '  v) 
E 


12 


'11 


-  I 


I 

12 


(5.5) 


(5.6) 


The  general  procedure  is  then  to  take  the  isotropic  solution. 


express  p  and  v  in  terms  of 


and 


by  means  of  (5.6)  and 


then  to  substitute  the  t  j  and  l  ^  (5.4).  The  plane-strain  solution 

for  relief  of  regional  stress  components  Oi  ,  03  ,  ^2  around  a  hole  of 
radius  a  in  an  isotropic  me  Uum  gives  a  radial  displacement  of 


u  1-  2 

^  [  o1  h  o2  +  2  (a1  -  o2)  cos2e  +  4r12sin26  ]  . 

(5.7) 


The  factor  containing  tne  elastic  constants  can  be  replaced  directly 
by  -t1  by  means  of  (5.5)  and  substitution  of  j  from  (5.4)  gives, 
for  tne  transversely  isotropic  medium 
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- 


a~=  "  (  bn  “  *  ^ai  1  °2  '  2(a1"a2)cos20  4  4T12sin29  1  .(5.8) 

3  3 

This  result  is  obtained  on  the  assumption  of  plane-strain,  in 
particular  that  ez  -  0.  Since  that  is  true  only  when  the  stress 
components  are  related  as  in  (5.2),  a  correction  term  must  be  added  to 
(5.8)  to  take  account  of  the  more  general  case.  It  was  noted  in 
Section  3  that  ez  is  independent  of  position,  so  that  we  may  substitute 
in  the  third  equation  of  (5.1)  the  values  of  the  stress  components  at 
infinity  to  determine  that  is  (since  the  effective  values  correspond 
to  removal  of  the  regional  stress), 

ez  ‘  -(b13Vb13°2  +  b33a3>'  (5-9) 


This  strain  is  due  to  the  fact  that  the  axial  stress,  q  ,  is  greater  than 

that  given  by  the  plane-strain  condition  (5.2)  by  a  quantity  az  (say), 

and  we  must  now  use  equations  (5. 1)  to  determine  c£  and  the  additional 

strains  e°  and  e°  which  it  causes . 
x  y 

By  putting  =  oy  =  0  in  (5. 1)  we  find  that 

o  0,0 

e  -  e  =  b,„o 


e 

z 


b33°z 


and  hence 


o 

e 

x 


and  so,  from  (5.9)  ,  in  terms  of  the  regional  stresses, 


o 

e 

x 


-b13CT3 


a2). 


(5.10) 


o  o 

Since  e  =  e  ,  the  additional  strain  is  independent  of  0  and  so 
induces  a  diametral  deformation  of  the  same  value  at  the  hole. 
Addition  of  (5.10)  to  (5.8)  leads  to  the  final  result  for  stress-relief: 
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(5.11) 


(5.12) 


(5.13) 


It  should  be  noted  that,  for  0  =  0,  the  x  and  y  directions  are 
defined  only  to  the  extent  that  they  are  normal  to  the  axis  of  elastic 
symme  try. 


5.2  Case  j  =  —  • 

As  noted  in  I,  the  plane-strain  problem  in  this  case  has  been 
discussed  by  Green  and  Zerna  (1954),  although  the  displacement 
solution  has  not  been  given.  However,  it  is  possible  to  use  the  complex 
potentials  given  by  them,  as  stated  in  1(5.31)  and  1(5.32),  with  the 
general  expression  lor  the  complex  displacement,  to  obtain  the  plane- 
strain  displacement  solution.  As  before,  the  stress-relief  solution  is 
completed  by  making  a  correction  for  the  non-zero  axial  strain. 

from  I  ( :  .3)  and  1(5.20),  the  stress-strain  relations  for  #  -  ~ 
are 
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ll 

b  a 

11  X 

b  a  -i 
13  y 

CM 

r—i 

X 

•<" 

II 

b,  _a 
13  x 

b  a 

33  y 

-t  b,  Oo 
13  z 

II 

WN 

b  „a 

12  x 

■  b  a 

13  y 

'  bi  ,a 

11  z 

2y  =  b  t 
yz  44  yz 


2y  -  2(b.  -b.  J  t 
zx  11  12  zx 


2y  -  b  t 
xy  44  xy 


expressed  in  the  coefficients  of  1(1 .3) . 


For  plane-strain,  e  -  y  =  y  =0 

z  yz  zx 


q  =  -  r-  (b. ,  a  b,  _o  )  , 
z  bu  12  x  13  y 


and  the  remaining  equations  become 


e  -  l,, a  l  0a  , 
x  11  x  12  y 


n  0  "1  ■tn  —  O  , 

12  x  22  y 


2y 


where 


xy 


"11 


b . .  t 
44  xy 


=  b 


11 


^12  b13 


11 

311 


bl2b13 

bll 


1  -  v 


i7 (1  -  V 
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1 

4 


-  'V  '-VV'T----  I.  UJ^ppp^ppp^|pp|«pp|p^!|p|pppp|ppp||| 


*\v 


33  b 


13 

11 


1/M  . 


(5.17) 


In  the  notation  of  Green  and  Zerna  the  complex  displacement 
due  to  this  kind  of  plane  strain  deformation  in  a  transversely  isotropic 
medium  can  be  expressed  thus: 

u  4  iv  =  (z j )  -  Pjfl’tZj)  +  6 2uu ’ (z2)  P2UU'(Z2^  '  (5.18) 

where 

61  -  v1)B2-<1-Yi)e1  . 

62  ~  ^  J‘  Y2^1~  ^  "  Y2^2  ' 

\  =  <!  +  Yi>%  ’  (1'V  S!  ’ 

\  =  f1  ’  y2>»1  "  Y2)  @2  ' 


Yj  a.  +  1  '  "  *12  "  *22  aJ  * 

J 

al  =  *12  +  ^b44  "  t  (  *12  ■*  ^b44J  "  ^2*22  ^  ' 
a2  =  *12  +  ^b44  "  l  ^12  H  ^b44*  "  *12*22  ^  ' 

and 

zj  *  2  +  v  -  1  =  u2  ■ 
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(5.19) 


i 

1  I 


1 


| 


rpyiy’S'V'  « >■  . 


For  problems  involving  circular  boundaries  the  z  are  more 
conveniently  expressed  in  terms  of  new  variables  C  and  C9  such 
that  1  1 

i 

and  then  ~'*2  =  ae*^  0n  boundary  (when  z  =  ae^).  The 
complex  potentials  for  the  present  problem  are  given  in  equations 
1(5.31)  to  1(5.33),  which  are  here  restated  in  the  form: 


Q'fZj)  -  f  (  Cj)  <  uu*(z2)  =  g  (  C2)  , 


where 


f(C)  =  Ac  -  [(l  -y1v2)a  ^  (i-  y2y2)A  ]  “cTv- r7T 

1  2 


9(0  -  A'  5  h  [  (1-YjVj)  A  +  Cl-  Vly2»  A- 1 


32A  = 


(  a  +  l)2(a9-i  l)2  ? 

2 - 2 -  L(1'y2)  °2)  '  2y2  (al~  ] 


2  2 
®1  “a2 


+  ^«V1,T12  ' 


32A' =  - 


(a  -i  l)2(cy  l)2 

2  2 
al  "°2 


[  (1  -•  Yj)  (a^  c^)  +2Y1(a1-  c2)  ] 


2i  .  ,  .2 

+  ^(a2+1)  T12  * 


On  the  boundary  ,  these  equations  lead  to  the  result 


(5.20) 


(5.21) 


(5.22) 


f I  (z  x )  =  f  (ae10) 


A  19 

Aae 


uu'(z2)  =  g(ae^ 


[(i-YjY2)a,  (i  -  y2y2)  A*  ] 


ae 


-ie 


Y1  "  Y2 


(5.23) 


..  ie  YjYj)  A  1  (1-  Y1Y2)  A'  _i0 

-  A  ae  i  - ae 

Y1  'Y2 


and  so,  from  (5. 18) , 


V  1U0  -i6  (u  +  iv) 

-  =  e  - 

a  a 


*  6,A  62A'  - 


3j[  (i-YlY2)  A't1"  Y2Y2^  A'^  '  P2[(1-Y1Y1)A-i(1-Y1Y2)A’] 


Y1  -  y2 


(5.24) 


-218  f  -  -  61ld-V1Y2)At,<l-Y2Y2)A,]-62[(l-v  v)A<(1-y  y  )A  ] 

■'  e  j  p  A+  p2A'  -  - 

l  Yi  -  Y2 

Upon  substitution  of  expressions  (5.23)  into  (5.24)  the  following 
result  is  obtained  a'fter  some  manipulation: 


u  +  iu„ 
r _ 0 


,2  j  11+(y1  +  Y2)Z-yJy2](01^H2(y1  +  y2)(^-02) 


(1  ” Y x )  (!-Y2) 


-2i(l-Y1Y2)(Y1+Y2)  t12 


+  e  2l^  2 (y1  YjHOj’’  c^>+2 (1  -i  Yj Y2)  (o^-Oj) 


(5.25) 
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Now  from  the  definition  of  and  Y2  in  (5.19),  they  are  either 
both  real,  or  complex  conjugates.  Tne  separation  of  (5.25)  into 
real  and  imaginary  parts  then  follows  easily,  and  we  have  for  the 
radial  displacement  due  to  plane-strain: 


-  A’  +  B‘  cos20  t  C1  sin2  0  , 


(5.26) 


where 


A'  =  - 


1  (V  Y2)2~  Y1  V2  1  (ai  *  °2>  *  2(y  v2)(0r  g2>! 


a  -  y/u  -  v2)2 


B'  -  - 


2t22l<VY2)(V  °2]i  (1'Y1Y2)(  V'V  ] 


(i  -  v/d  -  y2)2 


(5.27) 


C’  =  - 


4VI  '  Y1Y2>  t12 

U  -  Vj)2  0  -Y 2)2 


Briefer  expressions  may  be  obtained  (as  in  the  stress  solutions  of  I) 
by  using  the  constants 

2  2  t — — ■ 

1  1  1  /  “  “  \  1  f  n  /1  .1  \  •  ir- 


1  =  ,  k2  "  ^  (  °2  )-  k3  =  ^  2^ki  +  k2^  =  ^1  ^  *  ^5'28^ 


Tnen 


A'  =  -  H22  j  [  (ki~  1)  +  k;J(k1+  1)]  (o^  0^  (k1~l)(k1+  k^-i  1)  (ctj  -  o2)  , 

b'  -  -  h22(v  v  l)l  <V  1)(or  °2)’  <v  1)(ar  a2)/'  (s-29) 

C'"-*<22k2(kl‘k3J  1]  t12  ' 

By  substitution  of  values  equal  in  magnitude  but  opposite  in 
sign  to  the  regional  stress  components  in  the  third  of  (5. 14)  we  iind 
the  constant  axial  strain. 
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(5.30) 


e 

z 


=  -  <b12°l 


b13°2 


b 1 1 °3) 


and  calculate  the  additional  radial  displacement  caused  by  this 
departure  from  plane-strain.  Setting  ax  =  a  =  0  in  (5.14)  we 
find  the  relations  between  this  and  the  additional  strain  components 

°  ,  ° 
e  and  e  : 


o 

e 

x 


.  o 
=  b12°x 


o 


e 

y 


b  o° 
1 3  z 


e 

z 


(5.31) 


However,  since  we  are  interested  in  the  radial  strain  in  a  direction 
at  an  angle  0  to  the  x-axis  we  use  the  transformation 

2  2 

e  =  e  cos  0  •  e  sin  0=Me  ■>  e  -  (e-e  )cos2  0  ] 
rx  y  ^xyxy 

so  that,  from  (5.31), 


o  _  z 
£r  "  2b. 


11 


1  b12  b13  (b12'  b13>c°s2£lJ 


At  the  hole  this  corresponds  to  the  comparative  radial  displacement 

O  ,  . 

u^  /  a,  and  so,  using  (5.30), 


o 

u  =  - 
r 


IT-  <b12°l+  b13a2  +  bll°3)lV  V(b12-b13>  COs2e  1  • 


11 


(5.32) 


By  means  of  I (i .  3)  this  expression  can  be  put  into  the  form 

_V=  '4^  |(vl  +  v2)2(°l  '  (vl  '  v2>  (V°2)~2(vl  '  v2>  °3 


h  [(  Vj  -  ^MOj+c^)  -■  (v1-v2)2(o1-^)-2(v1-v2)  o3]cos2ej  (5.33) 
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which  can  then  be  combined  with  the  plane-strain  solution  given  by 
(5.26)  and  (5.29)  to  obtain  the  complete  solution: 

u 


■f  =  A90  +  B90COs2e  '  C90Sln2e  ' 


(5.34) 


where 


A90  = 


=  -*fVk>' 


(v,+  vj* 


»  '  k3<kl 


...i  •  } 


2  2 


(  V1  "  V2  )  Vl^ 

'[<22(kr1,(kl+k3"1)  (VC2 )ji  -217  ° 3  • 


B 


90 


•i(i22<kr1,(ki 


2  2 
V1  _  v2 

+k3+1)  "  (  V  °2> 


(5.35) 


r  (vv2}  i  vrv2 

['t22^kr1)(kl"k3+1)  "  Ej  j  (o1-CT2)  +  2Ej 


v,  -  vr 


3  ' 


C90  "^22k3(kl+ k3  +  1)  T12. 

The  constant  can  exPressed  in  terms  of  E1 ,  v1  and  k^  thus: 

2 


1  -  v, 

'22  ~  2 
L  F  k 

Er  i 


(5.36) 


5 . 3  Isotropy 

The  result  for  an  isotropic  material  may  be  obtained  by  letting 
E1=E2  =  E,v1  =  v2  =  v,  k1  =  1  and  k3  =  2  in  (5. 13)  or  (5.35). 

These  expressions  then  reduce  to 
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A 


B  ‘  -  |  (1  -  v2)  (  Oj  -  o2)  . 
C  -  -  |  (1  -  v2)  t12  . 


(5.37) 


which,  of  course,  are  the  coefficients  in  the  expression 

u 

—  =  A  +B  cos0  -f  C  sin0  .  (5.38) 

6 .  Application  of  the  radial  displacement  formulae 
6 . 1  Development  of  procedure 

In  each  of  the  expressions  for  the  radial  displacement,  (4.17), 
(5.11),  (5.34),  the  result  depends  on  only  three  parameters  which  are 
constant  for  any  particular  hole.  No  matter  how  many  measurements 
of  ur  are  made  at  various  directions  0  in  that  hole,  it  is  not  possible 
to  gain  any  more  information  than  that  of  the  values  of  the  parameters 
A0 ,  ,  C0  .  However,  in  order  to  determine  the  six  components  of 

regional  stress,  it  is  necessary  to  determine  the  values  of  six  parameters. 
To  do  so,  measurements  must  be  made  in  a  second  hole  with  a  different 
orientation.* 

In  each  hole,  a  minimum  of  three  diametral  strain  measurements 
is  adequate  but  more  may  be  made  to  improve  accuracy  (see  Panek,  1966, 
for  procedure  in  the  isotropic  case).  Denoting  the  three  measurements 
by  suffixes  1,2,3  we  have 

Ur  1 

— 1 —  =  A  +  B  cos  20  ,  +  C  sin  20 .  , 

a  1  1 

u  , 

— =  A  +  B  cos  20o  +  C  sin  20o  ,  (6.1) 

a  2  2 

u  „ 

-J—  =  A  +  B  cos  20o  +  C  sin  20»  . 

a  3  3 


*  In  some  circumstances  measurements  in  a  third  hole  may  be  necessary: 
see  note  at  end  of  Section  6.2. 
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The  determinant  of  the  coefficients  of  A,  B,  C  in  (6.1)  is 
4  sin  (Oj  -  )  sin  (S^  -  O3  )  sin  (63  -  Sj)  showing  that  these 

parameters  may  always  be  determined. 

When  measurements  in  two  holes,  with  orientations  0j 
ar.d  02  (say),  have  been  made,  six  parameters ,  Ap  Bp  Cp  A2 , 

B2 .  C2  ,  are  available  for  determination  of  the  regional  stress 
components.  For  definiteness,  we  suppose  that  0  0j  i  02  ^  n/2  , 
and  so  indicate  the  place  of  the  exceptional  values  of  0  in  the 
scheme . 

It  is  important  to  remember  that  the  regional  stress  components 
in  the  form  >  '^3'  T3l  are  referred  to  axes  fixed 

on  the  axis  of  the  hole  and  hence  will  change  value  as  the  orientation 
of  the  hole  is  changed.  Before  proceeding  to  solve  the  six  equations 
for  the  regional  stress  components  the  formulae  for  A^  ,  B^  ,  must 
be  expressed  in  terms  of  components  referred  to  a  single  system  of 
coordinates.  There  are  several  reasonable  systems,  but  it  seems 
preferable  to  choose  one  in  which  one  axis  is  an  axis  of  elastic  symmetry. 
We  denote  the  chosen  system  by  (  x^  ,  x^  ,  xR  )  and  the  components 
of  regional  stress  R...  Then  we  choose  Ox^  to  be  anaxis 

of  elastic  symmetry,  Jso  that  it  is  identical  with  Ox-<  (Fig.  3).  We 
could  choose  Ox^  to  be  a  geographical  direction,  or  to  be  identical 
with  the  Xj -direction  defined  for  either  of  the  holes  (but  not  by  both 
unless  both  holes  lie  in  a  plane  through  an  axis  of  elastic  symmetry) 

(Figs.  1  and  2) . 

Let  j  be  the  angle  made  by  the  Xj  -direction  of  the  first  hole 
to  Oxj  ,  and  \|r 2  the  angle  made  by  that  of  the  second  hole.  The 
direction  cosines ,  d^J  =  cos  (j^Oxfl)  (m  =  1,2),  of  the  x, -directions 
defined  for  the  first  hole  (m  =  1)  and  second  hole  (m  =  2)  with  respect  to 
the  (xR  )  coordinates  are 

,m  ,  ,m  ,m  n 

d, ,  =  cos  if  ,  d.  =  sin 7  ,  d,  =  0  , 

11  m  12  m  13 

,m  ,  ,m  ,m  ,r 

d  =  -  sin  i}'  cos0  ,  d  =  cos  vmcos0  .  ,  d  =  -  sin0  ,  (6.2) 

21  m  m22  m  : .  2  3  m 

,m  ,m  ,m 

d,  .  =  -sin  j  sin0  ,  o.n  =  cos  i  sin0  ,  d,,„  =  cos 0 
ol  mm  il  mmoo  m 

R 

If  we  choose  Ox^  to  be  identical  with  Oxj  for  m  =  1  (say),  then  =  0  , 
with  consequent  simplification  in  that  case. 


If  we  denote  by  (m  =  1  ,2)  the  regional  stress  components 
referred  to  the  xJ  -  coordinates  and  x^  -  coordinates  respectively, 
tnen  they  may  be  expressed  in  terms  of  the  R.j  (the  regional  stress 
components  referred  to  tne  xB  -  coo  inates)  by  means  of  the 
transformations 


m 


a‘“  =  d..  d  R,  (summed  over  k,  i  ) 
ij  ak  j l  kl 


(6.3) 


Written  out  in  the  original  notation  for  the  regional  stress  components 
this  gives  the  expressions 


m2  2 

a,  =  R,,cos  lr  +  R00sin  i  +  2R,„sin;|i  cos; 

1  11  m  22  fm  12  Ym  m 


m  2  2  22  2  2 

=  R,,sin  ill  cos  0  +  Rnocos  ill  cos  0  +  R_sin  0  -2Ri0sin;  cosli  cos  0 

2  11  m  m  22  Tm  m  3  3  m  12  m  m  1 


m 


-2R-„cosili  sin*  cos0  +2R,_simli  sin0  cos0  , 
23  m  m  m  13  m  m  m 


m 


2  2  2  2  2 
a‘“  =  R,,sin  tfi  sin  0  +  R00cos  i)i  sin  0  +  R,  .cos  0 

3  11  fm  m  22  ym  *m  3  3  *m 


-2R10sinili  cos\ii  sin  0  +  2R00cosil'  sin0  cos0  -2R,,.sini|i  sin0  cos0 

12  m  m  m  23  m  m  m  13  m  m  1 


m 


m  ,22. 

t,  =- R, , sinili  cos  ill  cos0  +  R„_simli  cosili  cos0  +R,„cos0  (cos  i 1  -  sin  li  ) 
12  1 1  m  m  m  22  m  m  m  12  m  m  m 


-  R.  simli  sin0  -  R_,, cosili  sin0 
23  Ym  m  13  m  m 


(6.4) 


m 

T23  " 


R,,sin%  sin0  +  R00cos  ill  sin0  cos0  -R._sin0  cos0 
1 1  m  m  iu  Z  m  m  m  j  o  m 


m 


2  2 

-  2R,0sinw  cosili  sin0  cos0  +  R0.cost|i  (cos  0  -  sin  0  ) 

12  m  m  m  m  23  m  m  m 


2  2 

+  R,  .sinili  (sin  0  -  cos  0  ), 

13  m  m  m 


m 


2  2 

t“‘  =  -  R, ,  sinili  cosili  sin0  +  R„  .sinili  cos  Or  sin0  +Riosin0  (cos  ill  -sin  1/  ) 
13  11  m  m  m  22  Tm  m  m  12  m  Ym  ni 


+  R„.  sin  it  cos0  +  R.  .cosili  cos0 
23  m  m  13  m  m 
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The  exceptional  cases  =  0,  ^  =  ^2  bring  simplifications, 
particularly  if  ft  or  ^  are  chosen  conveniently. 

(i)  =  0;  =  0  .  In  this  case  the  xj  -  direction  is  chosen 

arbitrarily  in  the  plane  of  symmetry  and  the  Xj  -direction  is 
chosen  to  be  identical  with  it.  Then,  from  (6.4) 


1  1  1 
a  =  R.  ,  ,  o  =  R00 ,  =  R,  „ 

1  11  2  22  3  ao 


T 1 2  R 1 2  '  T23  R23'  T 1  3  R1 3 ' 


(6.5) 


and  the  formulae  (5.13)  become 


Ao  ”  E1  ^11  +  R22  "  V2 R 3 3 ^ ' 


Bo  ~  E1  ^  ~  V2  ^  R1 1  R22^' 


Co  ~  “  El  ^  "  V2  *  R12  ' 


(6.6) 


(ii)  0  ,  =  /  ■  4r  =  0  .  This  choice  of  i|l  means  that  the 

Z  Z _ Z  2 

D 

x^  -  direction  is  in  the  plane  of  cross-section  of  the  hole, 
and  the  xB  -  direction  coincides  with  the  axis  of  the  hole. 
Since  the  x|  -  axis  for  0j  =  0  in  (i)  above  can  have  any 
direction  in  the  plane  of  elastic  symmetry,  the  choice  of 
|  =  0  in  each  case  does  not  prevent  the  selection  of  a 
common  xB  -  coordinate  system.  From  (6.4),  we  have  then 


Q1  =  Rll' 


°2  R33' 


”  R1 3 '  T23 


R22 ' 


T  1 3  R12  ' 


a'nd  the  formulae  (5.35)  become,  using  (5.36), 


(6.7) 


'-“**■'*“■'’****'*•*■±1  «;-»*«* 


'90  =  '«!  I  VLll(kr1)2  tk3(kr  1)1  +(V''2)2/<Rn+R33> 


2  2  )  V1^V2 

«J  —  (kr1,(kitki+1)+vi-^/(RirR33)+ ^  r22  • 


/  ‘i 


■f  (krl)(k1+Vl)+v^-v2 


*R11  +  R33^ 


(6.8) 


4E.  1  .2 


|  ^  <k. +  1)(k1+k3  +  1)+(vrV2jl<R.rR33) 


v.O  R 

2E  22 


U  -vt)k3  <k1*te3  *  1) 


'13  1 


Unfortunately,  we  cannot  make  appropriate  substitutions  (as  we 
have  for  0  =  0,  n  /  2  )  in  the  expressions  A^,  B  ,  (4.18),  for 

0  <  0  <tt/2  ,  to  obtain  formula  with  explicit  dependence  upon  the  R.., 
because  their  dependence  on  the  regional  stress  components  arises  in 
a  complicated  manner  through  the  a^Q  and  (and,  more  simply,  through 
ez  ).  The  dependence  of  the  a^  on  '  °2  ~12'  Tl  '  T23  can  *3e 

found  only  by  solving  equations  1(4.2),  and  that  is  not  practicable  until 
numerical  values  can  be  assigned  to  the  coefficients  (calculated  from 
the  elastic  properties' of  a  particular  material  with  the  appropriate  0  in 
each  case).  The  A^  can  then  be  calculated  from  1(4.12).  The  component 
c  is  introduced  through  ez,  which  is  given  in  terms  of  stress  components 
by  (3.2).  In  a  practical  case,  a  ,  Av  and  ez  would  then  be  expressed 
in  terms  of  the  ,  using  the  appropriate  values  for  the  angles  0m  and 
i|/m  for  each  hole  in  (6.4). 

The  elastic  parameters  G  and  also  depend  upon  0  ,  and 
the  values  appropriate  to  0^  ancl  0  must  be  calculated.  All  the  equations 
are  linear  in  stress  components  and  when  all  the  substitutions  have  been 
made  we  obtain  expressions  A j,  B^.Cj  for  0=0  ,  v=  ,  j.  linear 
and  homogeneous  in  the  stress  components  R^  ana  a  similar  set  A ^  , 

B  ,  C  for  0=0  ,  if)  =  ^  If  we  denote  byJ  A|  etc.  the  experimental 

vllues^determined  through  equations  of  the  form  of  (6.1),  we  can  write 
the  result  as  six  equations  thus: 
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A  ■**1 


I 


alRll  +a2R22  +  a3R33  +a4R23  +a5R31  +a6R12  =  A! 


blRll  +  b2R22  "'  b3R33  ■*  b4R23  +  b5R31  +  b6R  1 2  B1 


C1R11  +  °2R22  "  C3R33  +  C4R23  +  C5R31  +  °6R  1 2  C1 

aiRll  +  a2R22  +  a3R33  a4R23  +  d5R31  +  a6R12  =  A2 

blRll  '■  b2R22  '  b3R33  ''  b4R23  +  b5R31  I  b6Rl2  =  B2 

C 1R1 1  +  °2R22  +  C3R33  +  °4R33  +  °5R31  +  °6Rli  =  °2 


(6.9) 


So  long  as  the  matrix  of  coefficients  a^  etc.  in  (6.9)  is 
non-singular,  the  equations  may  be  solved  for  the  R^,.  For  0^  =  0, 
ij/j  =  0  expressions  (6.6)  are  substituted  for  the  first  three  expressions 
in  (6.9).  If  0 2  =  n/2  ,  =  0  expressions  (6.8)  are  substituted 

for  the  second  three  of  (6.9).  A  complete  solution  cannot  be  obtained 
by  having  0^  =  0  and  02  =  n/2  simultaneously,  since  R23  occurs  in 
neither  (6.6)  nor  (6.8)  and  the  matrix  of  coefficients  of  the  six 
equations  is  singular:  by  omitting  an  equation  a  solution  could  be 
obtained  for  the  other  five  components .  It  does  not  seem  practicable 
to  determine  whether  there  are  any  other  conditions  for  which  the 
matrix  is  non-singular.  However,  the  possibility  of  this  occurring 
must  not  be  overlooked  and,  indeed,  steps  should  be  taken  to  ensure 
that  the  matrix  is  "well-conditioned"  in  practical  applications.  This 
could  be  done  by  evaluating  its  determinant  at  a  selection  of  possible 
relative  orientations  of  the  holes ,  and  rejecting  those  giving  the 
smaller  values.  One  would  be  especially  suspicious  of  holes  at 
right-angles  to  each  other. 


6.2  Summary  of  method  of  application 


The  constants  k  are  evaluated  in  terms  of  the  b  (assumed 
known)  for  a  pafrS of  values  0  and  09  ,  through  .  6).  The 


constants  Ks<  \y.  Hv  .  Yv  .  6y.  Yyy.  6yy,  n#v,  L  ,  L2y,  Gy. 

H.  (  v  =  1,2,3)  are  then  determined  for  each  of  0,  ana  09  by 


\  V  J.  ,  L,  t  yjf  U1C  LHCU  UC  LCi  1II1I1CU  iWi  CUUli  ^  j  UllU  ^  O 

means  of  1(2.6'),  1(2.17),  1(2.18),  1(2.20),  1(4.6),  1(4.87, 
1(4.9),  (2.8)  and  (4.16). 


-5  I 


^*S#‘PT*»  :r-«  w-»i 


(2) 


(3) 


The  equations  1(4.2)  are  solved  for  a^Q,  a2Q,  a30'  a30  in 
terms  of  a.  ,  a  ,  t,0,  Ti3'  T23'  us*n9  t‘ie  coefficients 
appropriate  to  ea<?h  of  me  values  0  j  and  02  in  turn.  These  are 
substituted  in  1(4 . 12)  to  obtain  the  values  of  Av(v=  1,2,3) 
(using  the  appropriate  coefficients  for  each  of  0j  and  02  )  in 
terms  of  o1 ,  o^,  t^,  T93.  The  values  of  the  krs  for  0^ 

and  02  are  used  in  turn  in  (*1.2)  to  obtain  the  values  of  ez  in 
each  case ,  in  terms  of  •  °2  •  °j '  T23  • 

R  R 

The  orientation  of  the  Xj  ,  x2~  axes  in  the  plane  of  elastic 
symmetry  is  selected.  The  xj  -  direction  is  defined  (Figs.  1 
and  2)  by  the  intersection  of  a  cross-section  plane  of  the  first 
hole  with  a  plane  of  elastic  symmetry,  and  i|r,  is  defined  as  the 
angle  it  makes  to  the  x^-  direction:  ^  is  denned  similarly. 


(4)  The  values  of  the  a^o,  a  ,  Ay  and  ez  for  each  of  0  j  and  02 
are  expressed  in  terms  of  the  by  means  of  (6.4),  using  the 
t|im  appropriate  for  each  hole . 


(5)  The  resulting  expressions  for  a  ,  aVQ,  A  ,  Av  and  ez  for 
0  =  0^ ,  ijf  =  \ji^  are  substituted  in  (4.18)  with  coefficients 
appropriate  to  01  and  form  the  first  three  expressions  in  (6.9). 
A  similar  process  for  0  =  02,  i|f  =  i|i2  gives  the  second  three 
expressions.  If  0j  =0,  the  expressions  (6.6)  can  form  the 
first  three  of  (6.7);  if  02  =  n/2'  (6.8)  can  form  the  second  three 
of  (6.9),  so  long  as  the  xjj*  -  coordinates  are  appropriately 
selected. 


(6)  The  determinant  of  the  coefficients  of  the  R^  in  the  resulting 
set  of  expressions  is  evaluated  to  ensure  chat  it  is  not  zero 
or  very  small. 

(7)  If  this  test  gives  a  satisfactory  result,  stress-relief  tests 

are  made  in  boreholes  at  the  selected  angles.  Diametral 
measurements  are  made  in  three  directions  to  obtain  a  set 
of  equations  of  the  form  (6,1)  for  each  borehole.  From  the 
first  hole  we  calculate  dimensionless  parameters  A^ ,  B|  ,  , 

and  from  the  second  parameters  A| ,  b|,  C|  •  These  form  the 
right-hand  sides  of  equations  (6.9)  which  can  now  be  solved 
for  the  regional  stress  components  R. . . 

Note  If  the  ground  is  pnly  weakly  anisotropic  it  is  probable  that 
the  equations  (6.9)  may  not  be  well-conditioned  since,  as 
shown  by  Gray  and  Toews  (1967),  and  also  by  Bonnechere  and 
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Fairhurst  (1968)  with  reference  to  the  "doorstopper"  method, 
the  six  equations  obtained  from  any  pair  of  holes  are  not 
independent  wnen  the  ground  is  isotropic  and  measurements 
in  three  holes  are  necessary.  The  arguments  of  the  above 
authors  do  not  apply  to  ground  which  is  only  transversely 
isotropic,  although  instances  of  dependence  between  the  six 
equations  (6.9)  other  than  that  noted  in  Section  6.1  (0^  =  0, 

02  =  n/^  may  exist.  If  the  matrix  of  coefficients  of  (6.9)  is 
zero,  or  close  to  zero  (either  because  of  weakness  of  the 
anisotropy  or  because  of  limitations  in  choice  of  hole-direction), 
it  will  be  necessary  to  use  a  third  hole  (0  =  0^  ,  i|i  =  ^  )  in 
addition,  and  the  results  may  be  processed  by  methods  discussed 
by  Panek  (1966)  and  Gray  and  Toews  (1967). 

7.  Strain  in  the  wall  of  a  stress-relieved  borehole 


Development  of  very  small  strain-gauges  makes  it  feasible  to 
measure  the  strain  in  the  wall  of  a  borehole  as  the  stress  is  relieved. 
This  technique  could  be  used  to  augment  information  obtained  from  the 
older  method,  and  possibly  to  allow  determination  of  the  regional  stress 
tensor  from  measurements  in  a  single  borehole. 


The  components  of  strain  in  the  wall  of  the  hole ,  r 


(r,  9,  z)  coordinate  system  are  ea  ,  e  ,  y 


constant  and  is  given  by  (3 


uic 

.2),  wh 


ile 


bz 


Of  these 


2  2 
en  =  e  sin  9  -  2  y  sin9cos9  +  e  cos  9  , 
9  x  xy  y 


Y„  =  v  cos6  -  y  sin9 
9z  yz  xz 


7.1  General  case,  o  <  0  <  n/, 


in  the 
is 


(7.1) 

(7.2) 


In  this  case  e  ,  e  ,  y  ,  y  ,  yxz  are  given  by  equations 
(2.3)  to  (2.8)  and  €  byy(3.Zj,  with  the  W..(z  )  =  -  [W  (z  )  ] 

as  given  by  (4 . 9) .  From  (7.1),  (2.3),  (2 .4)  and  (2 . 5) ,  we  find  that 
3 

en=S  [  M  W  (z  )  +  M  W  (z  )  ]  +  (  K.  sin20  +  Kocos20  )  e  ,(7.3) 
9—1  v  v  v  v  v  v  1  2  z 

v  =  1 

where  2 

M  =  L.  sin  9  +  L„  cos  9  -  L,  sin9cosB  ,  (7.4) 

v  lv  2v  6v 
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L  ,  L„  and  L  being  given  by  (2.8).  For  the  boundary  of  the 
hoYe,  rv=  a,  we  Yise  the  value  of  [  Wv(z^)  ]  ,  given  by  1(4.13), 

so  that  then 


W  (z  )  = 
v  v 


(a  y  -  A  )  sin0  -  i  (a  y  +  A  )  cos  9 
vo  v  v _ vo  v  v 

sin9  -  X  cos0 
v 


(7.5) 


From  (7.2),  (2.6)  and  (2.7),  we  find  that 


Y  =  )  [NW(z)+NW  (  z  )  ]  +  K.  cosy  e  ,  (7.6) 

0z  -  .  v  v  v  v  v  v  4  z 

v  =  1 


where 


N  =  \  (L.  cosS  -  sin9  )  , 

v  4v  5v 


(7.7) 


and  L,.  being  given  by  (2.8)  ,  and  W^(zJ  by  (7.5),  for  r  =  a, 
In^>oth  (7.3 i'j  and  (7.6),  the  value  of  e  is  given  by  (3.2): 


ez  =  -(k130l+k23°2+k33C3  +  k34T23  >  ' 


(7.8) 


Equations  (7.3)  to  (7.8)  give  the  distribution  of  strain  in  the 
wall  of  the  hole  as  a  function  of  the  elastic  constants,  hole  orientation 
and  regional  stress  tensor.  The  dependence  of  a^Q  and  A  on  the 
components  oi  regional  stress  must,  in  practice,  be  calcuYated  by  the 
methods  indicated  in  Section  6. 

7 . 2  Case  0=0. 

In  the  case  0=0,  ^ / 2  >  the  plane  deformations  and  the  axial 
shear  ocfo. mations  arc  "uncoupled"  and  so  eg  and  Ygz  may  be 
calculated  as  separate  problems. 

Transformation  of  the  plane-strain  relations  (5.1)  to  the  (r,9,z) 
system  of  coordinates  gives  the  relation 

e.  =  b.  .0,  +  b  0a  b  ,0 
0  lie  12r  13  z 

but  since  c  =  0  on  r  =  a,  we  have 
r 


eeibaVbi3Von,=a 


(7.9) 


f, 


where  b  and  b  are  given  by  1(1 . 3) . 

i.  1  X  O 

The  values  of  and  oz  may  be  obtained  from  1(5.5)  and 
1(5.7)  with  the  signs  changed  (since  here  the  regional  stress  is 
being  removed): 

o  =  -  ^1  +  a2^  +  2  (°i  "  a2  )  cos29  +  4  Tl2sin28  ,  (7.10) 


a  =  -  cl  -  - —  [  2 (a.  -  a_)  cos20  +  4  t.  osin20  ]  .  (7.11) 

2  3  b33  1  2  12 

Substitution  of  (7.10)  and  (7.11)  into  (7.9)  give  the  result 
for  r  =  a , 


e0  =  -bil(ol+c^  "bl3a3  +2(bll  --r;)Ua1  -a2)cos29  +2T12sin20 

J  J 

(7.12) 

From  1(1.3)  this  may  be  written 


El£0  =  ‘^°1+  a2^+  V2a3  +  "  E  V2  ^  ^  (ar  a2^cos2e  +  2  Tl2sin29  1 

i 

(7.13) 

By  changing  the  sign  of  1(5.19)  (because  the  regional  stress 
is  removed)  we  have  the  component  of  shear  stress  in  the  wall  of 
the  hole: 


=  -2  (t„ „cos9  -  t.  osin0  )  . 
9z -  2 j  13 


(7.14) 


(From  (5.1)  (suitably  transformed)  we  have  the  relation 


2y8z  b44T0z 


and  so ,  from  (7.14) 


Yte“b44(T13sin9  -t23COs9)  ' 


or ,  using  1(1 . 3) , 


My„  =  t  _sin9  -  TOOCOS0. 
0z  13  23 


(7.15) 


(7.16) 
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is  given  by  equation  (5.9),  which, 


=3.  (7.17) 

Equations  (7.13),  (7.16)  and  (7.17)  give  the  components  of 
strain  in  the  wall  of  the  hole  directly  in  terms  of  the  elastic  constants 
and  the  regional  stress  components  referred  to  axes  fixed  on  the  hole. 

7.3  Case  0  =  71 r/ ' 

A  slightly  different  procedure  is  used  for  0  =  n/  2  in  orc*er  t0 
avoid  the  labour  of  transforming  the  equations  (5.14),  which  apply 
to  this  orientation  of  the  hole . 

The  stress  components  ox  ,  a  ,  Txy  are  derived  from  components 
referred  to  the  (r,  9  ,  z)  coordinate/ by  the  transformation 

2  2 

0  =  0  cos  0  +  o^sin  9  -  2t  .sinftcosO  , 
x  r  9  r0 


2  2 

o  =0  sin  9  +  a.cos  9  +  2  t  .sin9cos0  , 
y  r  0  r0 


2  2 

t  =  (a  -  aj  sin9cos0  t  .(cos  0  -  sin  9). 
xy  r  u  r0 


(7.18) 


The  remaining  component,  ez 
by  use  of  1(1 .3)  can  be  written: 


Elez=  V°1  +  °2>  -  t, 


In  the  present  case,  for  r  =  a,  or  ~  Trg  =  0  ,  and  Oq  is  given  by 
1(5.40)  and  1(5.41),  and  oz  by  1(5.46)  both  with  the  signs  of  the 
regional  stress  components  changed  (since  they  are  being  removed, 
not  applied) . 


From  (5.14) 


e 

x 


b,  a 
11  x 


bi  qct 

13  y 


Vz 


e 

y 


b.  .a 
13  x 


b  a 
33  y 


b.  ,a 
13  z 


b 


.  .T 

44  xy 


(7.19) 
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Substitution  of  (7.18)  in  (7.19)  and  use  of  (7 .  1)  gives  the  result 


4  2  2  4  2  2 

€  „  =  a  (bj  jSin  5  +  2b^  3sin  9cos  9  +  b^cos  9  +  b^sin  0  cos  9) 


(7.20) 


2  2 

h  o  (b,  „sin  9  +  b._cos  9). 

Z  1 2  13 


From  1(1.3)  and  1(5. 4b),  this  can  be  written 


F  F 

4  2  2  1  4  1  2  2 

E ,  =  a J  sin  0  -  2v„sin  Geos  9  +  —  cos  9  +  ~  sin  Geos  9 

109  2  M 


2  2  2  2  2 
(VjSin  9  +  v2cos  9  )  ]  +  fysin  9  +  v2cos  9)(o3  -  v  a  -  ) 


(7.21) 


where,  from  1(5.40  )  and  I  (5.41) 


H-aQ  =  -  (a^+  a2)[  k^k^  1)  +  (kJ-  1)  -  (k^  k3+  l)(k1~  1  )  cos  2  9  ] 

-(Oj-  (J2)(k  +  k3  +  1)  [  k  -  1  -(k  +  l)cos29  ]  +2  T12k3^i+k3+  1)sin29 


(7.22) 

H  =  (kj  +  2k2  +  l)-2(k^-l)cos20  +  (k^-2k2  +  l)cos229  (7.23) 

To  calculate  the  component  Ygz  we  use  (7.2)  in  conjunction 
with  the  relations  from  (5.14): 


V  =  £b  .  „t 
yz  44  yz 


yxz=(bn-b12>  Txz' 


(7.24) 


T 

XZ 


with  signs  changed. 


and  the  expressions  1(5.  bb)  for  t  and 
The  result  is  yZ 


uz 


(i  +  t) ( t  .  sin0  -  tToocos0) 

_ 1  J  Z  o _ 

2  2  2 
sin  u  +  t  cos  9 


[  ib44cos29  -■  (t1!  x -bi2)sin2S  ]  .  (7.25) 


where  ,  from  1(5 . 08) 


2 

t  = 


44 


,(bn-bi2> 


2M(1  +v  ) 


(7.26) 


Using  1(1.3)  and  (7.26),  this  maybe  written 


ElY0z  =  ^  +  V  ^  +  fc)  (  T13sin0  "  tT23COs9^ 


(7.27) 


The  axial  strain  is  given  simply  by  1(5.3),  which  may  be  written  , 
by  means  of  1(1 . 3): 

ElCz  =  -°3+Vl°l+V2  •  <7-28> 

Equations  (7.21),  (7.22),  (7.23),  (7.27)  and  (7.28)  give  the 
components  of  strain  in  the  wall  of  the  hole  (for  0  =  n/2  )  directly 
in  terms  of  the  elastic  constants  referred  to  axes  fixed  on  the  hole. 

7 . 3  Isotropic  Material 

The  results  for  an  isotropic  material  may  be  obtained  simply 
by  putting  E^  =  E2  =  E,  =  V2  =  v  ,  M=pi  =  £E/(1  +  v  ),  in  (7. 13)  , 
(7.16)  and  (7.17).  The  results  are 

2 

EeQ  =  -  (c^  +  c^)  +  vo3+2(1-v  )[  a2)cos29  +  2  T^sh^e  ]  ,  (7.29) 
EyQz  =  2(1  +  v  )  (  T13sin0  -  t23cos0  )  ,  (7.30) 

Ee  =  v  (  a  +  a_  )  -  a.  . 
z  12  3 
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(7.31) 


7.4  Measurement 


In  order  to  compare  the  results  of  this  section  with  measured 

values  we  require  the  normal  components  of  the  strains  calculated 

above  resolved  in  any  direction  tangential  to  the  wall  of  the  hole. 

Let  be  the  normal  component  of  strain  in  the  wall  in  a  direction 

at  an  angle  uu  to  the  positive  9 -direction  at  any  point  then  its 

value  in  terms  of  e  ,  e  ,  y  at  that  point  is  given  by  the  transformation 

0  z  9  z 

2  2 

e  =  e„cos  uj  +  e  sin  uu  +  2  y.  sinuu  cos  uu  .  (7 . 32) 

ju  9  z  0z 


If  a  triple  straii  -gauge  rosette  were  used,  ^with  gauges  at^ 
equal  angular  intervals ,  then  the  three  angles  0,  /„  ,  /„ ,  or  A  , 

tt  Bn  Job 

/„,  /.  would  be  obvious  choices  for  uu  . 

L  b 
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